Rescattering for extended atomic systems 
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Laser-driven rescattering of electrons is the basis of many strong-field phenomena in atoms and 
molecules. Here, we will show how this mechanism operates in extended atomic systems, giving rise 
to effective energy absorption. Rescattering from extended systems can also lead to energy loss, 
which in its extreme form results in non-linear photo-association. Intense-laser interaction with 
atomic clusters is discussed as an example. We explain fast electron emission, seen in experimental 
and numerically obtained spectra, by rescattering of electrons at the highly charged cluster. 
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Laser-driven rescattering of electrons [1] is at the heart of 
strong-field atomic physics. The basic principle is very 
simple: a bound electron is released from an atom (or 
a negative ion [2]) with the help of the strong electric 
field of a laser, by which it is subsequently driven back 
to the ion. On return to the ion the electron may re- 
combine (emitting harmonic radiation [3] with the access 
to attosecond laser pulses [4]), ionize other electrons (in- 
ducing multiple ionization [5]) or may be backscattered 
(gaining high kinetic energy termed above-threshold ion- 
ization [6]). Due to the strong dependence of the tunnel 
probability on the field the release time is restricted to 
phases of the laser period with maximal electric field. In 
the following we concentrate on linearly polarized light 
which exhibits the most pronounced rescattering effects. 

Most direct evidence for the rescattering mechanism 
comes from measuring kinetic energies of the released 
electrons [7-10]. The momentum p an electron can ac- 
quire in an oscillating field /(t) = F cos{Lot) of strength 
F and frequency lo depends on the phase Lp^ = Luto at the 
time to of its release. At this time the vector potential of 
the field is Aq :— Asinipo = (F/uj) sin(po and this is ex- 
actly the momentum acquired, p ~ Aq. (We use atomic 
units.) Obviously, the maximum pmax — A = F/uj occurs 
at = 7r/2 and results in a kinetic energy iS'max = SE'pond 
with the ponderomotive energy E'pond := F"^ /Auj'^. How- 
ever, electrons with such high energies are rare [8] since 
they must be released from the atom when the electric 
field vanishes, which is very unlikely (see above). Elec- 
trons with energies even beyond 2i?pond are indicative for 
rescattering of electrons previously released. 

For atoms the rescattering process can be understood 
in classical terms by an electron elastically [28] scattered 
at a zero-range potential in the presence of an electro- 
magnetic field [11]; neglecting the Coulomb tail turned 
out to be of minor importance. In this case, and more 
generally, if the scattering time is much shorter than the 
laser period, the absorbed energy AE' is given by [12] 



l^E = 2Ap • A{t') 



(1) 



with Ap the momentum change at the time t' of scat- 
tering. The absorption is particularly effective for a 
large momentum change Ap, which is achieved through 



backscattering, and for scattering events taking place at 
a maximal vector potential (minimal electric field). 

For an extended scattering potential we will see 
that optimum energy absorption upon rescattering re- 
quires completely different conditions, namely forward- 
scattering at maximal electric field. 

This becomes obvious by analyzing the simplest situa- 
tion of scattering from an extended potential. It can be 
realized with an electron passing under the influence of 
an oscillating electric field centrally a spherical well. The 
one-dimensional problem is sketched in Fig. 1. Assuming 
a constant potential — V" < inside the extended scat- 
terer naturally generalizes the zero-range potential used 
in the context of rescattering at atoms [11]. Note, that 
similar to the discussion before Eq. (1) for the atoms, the 
time of crossing the potential boundaries may assumed 
to be short compared to the laser period. Therefore, the 
boundaries are idealized here as steps in the potential. 
This allows for the corresponding electron dynamics to 
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FIG. 1: Schematic picture of the scattering process. The 
upper part shows the electron's momentum, the lower part 
the potential and the value of the electric field (assumed to 
be homogeneous in space) as a function of time, respectively. 
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be solved analytically. The electron momentum 

Pa=Pa + {F/UJ) COs{ujt) (2) 

contains a constant drift Pa (shown by dashed lines in 
Fig. 1) and an oscillating part due to the laser field 
(shown by thick solid lines in Fig. 1). In order to conserve 
energy across the boundaries of the potential reached 
at times and U, respectively, the electron momentum 
must jump at those times. 

The motion is determined by the laser (with field 
strength F and frequency w), the potential (with width 
L and depth V), and the initial drift momentum pi, cf. 
Fig. 1. A given entrance time tj (or phase ip^ := Loti) fixes 
the exit time ti and phase ipi := ajt{. Then the final 
momentum, cf. Eq. (2), reads 



Pi = 



1/2 



(3) 



whereby the phases (fi and ipi are connected by 



Liv = (-A + VWT^) (<Pf - ^i) - (4 - ^0 • (4) 

In addition to Aa — Asini^o,, we have defined A'^ := 
Acosipa- The "transit" momentum p := V2V charac- 
terizes the depth of the potential. We assume Pi > A 
to guarantee monotonic motion of the electron and avoid 
multiple passing of the potential border. 

The final momentum pf cannot be expressed explicitely 
in terms of (pi since Eq. (4) is essentially non-algebraic. 
However, it is instructive to optimize pf with respect to 
(fi and ipf. Then, Eq. (4) reveals the optimum width 
of the potential for extremal energy absorption under 
given laser light. One obtains extrema for ipi^f = ±7r/2 
mod 2tt and the maximum pf for ipi — —tt/2 and ipf = 
-|-7r/2 with the electron passing the center of the potential 
at times between —tt/oj and +7r/w. Electron momentum 
and potential width read 



Pf,ii 



= -A 



2A + ^{pi-Af+p^^ -p2 



A + ^^{pi-Ar+p^ 



1/2 

(5) 

(6) 



In the limit of a deep potential, i.e., p^Pi and p^ A, 
both expressions simplify to 



Pf ,inax 

L 



I^Ap 



max ^ P 

00 



(7) 
(8) 



Obviously the optimal potential width L in Eq. (8) corre- 
sponds to the distance an electron travels with momen- 
tum p during half a cycle tt/o; of the laser pulse [13]. 

At optimal width Lmax the maximum electron momen- 
tum achievable depends on both, the laser field amplitude 
A and the depth of the potential (through p in Eq. (7)). 



This is different from above-threshold ionization of single 
atoms where the maximal electron momentum is deter- 
mined by the laser field only, namely Pmax — A (or 
i^kin,max = lOi^pond) [6]. It should be mentioned that 
electrons with such high energies are typically a few or- 
ders of magnitude less abundant than low-energy elec- 
trons [9, 14]. 

The energy absorption is maximized by an increased 
momentum during one half cycle of the pulse. In a full 
laser cycle a free electron absorbs as much energy as it 
looses. If, however, the drift momentum is increased by 
p for just one half of the cycle there is a net energy ab- 
sorption. Using Eq. (2) it is given by 



-I 



AE = / f{t)P{t) dt = Fp 



+tt/2u 

I 

-7r/2w 



cos(a;i)di = 2Ap, (9) 



in accordance with the result in Eq. (7). Whereas the 
absolute change of the momentum in the field, being 2A, 
is independent of the drift momentum p the change of the 
kinetic energy Eq. (9) is proportional to it. Rewritten in 
terms of energies, we obtain for the maximum energy gain 
through rescattering from extended systems the central 
result of this paper: 



ae = 4,/e^Vv. 



(10) 



The acceleration of the electrons depends on both, the 
laser field strength A and the depth V of the scattering 
potential. Interestingly, this situation is akin to the so- 
called powered swing-by (or gravity-assisted maneuver) 
of spacecrafts [15], There, thrust for accelerating, decel- 
erating or redirecting the spacecraft is applied only for 
short intervals of time and not in an oscillatory fashion 
like in the case of a laser. However, similar to the situ- 
ation considered here, thrust is most effectively applied 
when the space craft has high momentum, which is the 
case at the perihelion of a swing- by at a planet. 

As mentioned before, for a given potential the phases 
(fi and iff are linked through Eq. (4). For the case of 
deep potentials, when p^pi and p^A, this equation 
simplifies to Lui fa p{(pi — ipi) . This fixes the difference of 
the phases and one may write the momentum explicitely 

1/2 

aspf((/5) = 2 [Ap cos (<^/2) sin (La;/2p)] ' with(fi ■— (pf+ 
ipi the only parameter left for optimizing pf . For (p = 0, 
i. e. when and tpi are symmetric with respect to the 
field maximum, the momentum reads 



Pf = Pf .max 



sm 



2Ln 



-,1/2 



(11) 



with the maximal momentum from Eq. (7) and the opti- 
mal system width Lmax from Eq. (8). 

To illustrate the relevance of the rescattering mech- 
anism in extended systems we present in the following 
results of microscopic calculations for rare-gas cluster ex- 
posed to intense laser radiation [16]. The theoretical ap- 
proach describes the laser-driven electronic nano-plasma 
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FIG. 2: (Color online) Kinetic energy distribution of electrons 
for two laser pulses of 400 fs duration and peak intensities 
of 1x10^^ W/cm^ (a) and 5x10^^ W/cm^ (b), respectively. 
Shown are the results for three different cluster sizes Xeias, 
Xeiisi and Xegogs, respectively. The first two are shifted up- 
wards for better visibility. Each set of points is fitted (solid 
lines) by an exponential curve exp(— _E/i5kin). Note that the 
energy scale is given by the corresponding ponderomotive en- 
ergy -Epond which is proportional to the laser intensity I. 



and the ionic explosion dynamics by means of classical 
equations of motion. It has been successfully applied to 
study, e.g., absorption mechanisms for a wide range of 
clusters sizes and laser pulses by various groups [17-20]. 
Figure 2 shows kinetic energy distribution of electrons 
as obtained for xenon clusters of different sizes and var- 
ious laser pulse parameters. The calculations are rather 
expensive since the the electrons have to be propagated 
for a long time (typically a few picoseconds) in order 
to obtain converged results for their final kinetic energy, 
which otherwise would be spoiled by the large, long-range 
and time-dependent attractive potential of the exploding 
cluster. All distributions show an exponential behavior in 
accordance with observations from experiments for some- 
what larger xenon clusters [21, 22] as well as silver cluster 
of similar size [23] . Fitting an exponential shape to these 
distributions yields parameters i?kin which are listed in 
Tab. I. They strongly depend on the cluster size and the 
laser pulse, but reveal a clear trend: The larger the clus- 
ter the faster the emitted electrons. This trend originates 
in the deeper potential for larger clusters, i. e. larger V in 
Eq. (10), which leads to a stronger acceleration at rescat- 
tering. 

Note, that the energy can exceed the ponderomotive 
energy considerably as seen in particular for the largest 
cluster Xegogs. Knowing the cluster potential from the 
simulation we can estimate the electron energies with the 
rescattering model of Eq. (10) and list them for com- 
parison in Tab. I. They agree surprisingly well consid- 
ering the simplicity of the model and the fact that the 
electron spectrum of the microscopic calculations (Fig. 2) 
contains also all electrons released directly. We attribute 
the agreement to the fact that the exponential tail is 
due to the fast electrons which are dominantly emitted 
by rescattering. Additionally, most of the electrons are 
ejected at the resonance of the cluster [22, 23] where the 
acceleration is optimal. This is certainly not the case 



TABLE I: Electron energies in units of the ponderomotive 
energy -Epond from xenon clusters of three sizes and two dif- 
ferent laser pulses. The exponential fit parameter Skin from 
the microscopic calculations shown in Fig. 2 are compared to 
the value A_B of the rescattering model in Eq. (10). 
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cluster size A'^ 
laser pulse 
microscopic result 

£kin/-Bpond 0.78 0.27 1.98 0.60 5.43 1.74 
rescattering model 

AE/Epond 0.19 0.09 1.53 0.65 3.68 1.58 



for the smallest cluster, Xei35, considered here. It is al- 
most completely disintegrated at the time when the laser 
pulse reaches its peak. What we assume to determine 
the parameters for rescattering is therefore not valid and 
consequently poor is the quantitative prediction of the 
electron spectrum by rescattering for this cluster. 

Similar to the calculated electron energy distributions 
experimental spectra can be characterized by constants 
^-kin quantifying the exponential decay. They are shown 
for measurements [23, 24] of various cluster along with 
the theoretical results discussed above in Fig. 3. Clearly, 
they are larger than the ponderomotive energy, i.e., 
above the dashed line in Fig. 3. Note that corresponding 
data for atoms [9] are below this line. Even in cases where 
one observes a plateau energetic electrons from atoms are 
much less likely than for extended systems such as clus- 
ters. 

In contrast to the microscopic calculation we generally 
neither know the charge nor the radius of the cluster at 
the peak of the laser pulse. However, there are two ex- 
ceptions: (i) the pulse is very short [24] or (ii) the delay 
of a dual pulse is adjusted to induce resonant ionization 
[23]. In both cases one can roughly estimate the un- 
known cluster parameters which determine width L and 
depth V of the extended potential for rescattering. In 
case (i) one can neglect the cluster expansion, the radius 
of the scattering potential is the initial cluster radius R. 
The charge Q can be estimated from a simple over-the- 
barrier model Q — FE? [25]. Hence, with the potential 
depth 3Q/2R of a homogeneously charged sphere we 
can by means of Eq. (10) determine the electron energies 
(open diamonds in Fig. 3) in reasonable agreement with 
the experimental values. In particular the ratio between 
the larger and the smaller cluster size is well reproduced. 
Note that the experimental signal is a sum over clusters 
of different size and the laser focus. This may be crucial 
[26] and agreement with a simple model on the abso- 
lute scale cannot be expected. For case (ii) we assume 
that the cluster expands homogeneously over the delay 
between the (short) double pulses. Thus one has at res- 
onance Q/R^ = lli'^ [20]. The charge Q can be roughly 
assessed from measured final ion charges [27]. Assuming 
an average ion charge of two we estimate a value (open 



4 



1 1 IIMI| 
■ 


1 1 IIIM| 1 1 1 IMM| 


1 1 IMMj 

♦ 


: 


• 

O 

8 s 


o - 


~ 1 


1 1^ 1 





0.01 0.1 1 10 
ponderomotive energy [keV] 

FIG. 3: Electron energies as a function of the ponderomo- 
tive energy iJpond from various clusters. The filled symbols 
show Skin for experiments (red) and microscopic calculations 
(blue): Agiooo (square, [23]), Ariroo and Arasooo (diamonds, 
[24]), Xensi and Xegogs (circles. Fig. 2). The corresponding 
estimates AE from the rescattering model (10) are shown by 
open symbols. 

square in Fig. 3) slightly above the experimental one. 

The experimental and numerical examples of rare-gas 
clusters demonstrate that the simple rescattering mech- 
anism for extended potentials we have introduced can 
provide considerable insight into complicated many-body 



dynamics as it occurs in these clusters including semi- 
quantitative predictions. However, the analytical results 
from Eqs. (3) and (4) are far more general than the exam- 
ple of rare gas clusters may suggest and should, e. g., also 
describe non- linear absorption of laser energy in quantum 
dots. Moreover, rescattering in extended systems does 
not only provide conditions for optimum energy absorp- 
tion from the light, it also answers the opposite question: 
Given a certain rescattering potential and laser pulse, 
what is the maximum velocity of a particle which can be 
brought to rest (remains sticking in the scattering sys- 
tem) under the combined action of laser and potential? 
This situation is also described by Eq. (5) if one inter- 
changes the indices "f" and "i" and sets pi = 0. The 
result is given by Eq. (7), but now pi,max ~ "i^/Ap stands 
for the maximum initial momentum an electron can have 
and still will stick to the extended potential being most 
efficiently decelerated by the potential and the light. This 
constitutes non-linear photo-association in extended sys- 
tems. 

We acknowledge support from the Kavli Institute for 
Theoretical Physics (KITP) in Santa Barbara, where this 
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